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CONNECTION BETWEEN THE LATTICE BOLTZMANN EQUATION AND THE BEAM 

SCHEME 


KUN XU* AND LI-SHI LUO t 

Abstract. In this paper we analyze and compare the lattice Boltzmann equation with the beam scheme 
in details. We notice the similarity and differences between the lattice Boltzmann equation and the beam 
scheme. We show that the accuracy of the lattice Boltzmann equation is indeed second order in space. We 
discuss the advantages and limitations of lattice Boltzmann equation and the beam scheme. Based on our 
analysis, we propose an improved multi-dimensional beam scheme. 

Key words, kinetic theory, hydrodynamics, lattice Boltzmann method, beam scheme, Euler and Navier- 
Stokes equations. 

Subject classification. Fluid Mechanics 

1. Introduction. The method of lattice Boltzmann equation (LBE) [1, 2, 3, 4, 5] is a gas kinetics 
based method invented to solve mainly hydrodynamic systems described by the Navier-Stokes equations. 
The lattice Boltzmann equation is fully discrete in time and phase space. It is a drastically simplified version 
of the continuous Boltzmann equation [6, 7]. The lattice Boltzmann method has a number of computational 
advantages: simplicity of programming, intrinsic parallelism of the algorithm and data structure, and con- 
sistency of thermodynamics. These advantages of the lattice Boltzmann method are due to the following 
facts: first of all, the Boltzmann equation has a linear convective term; second, the velocity space is re- 
duced to a set of very small number of discrete velocities in the lattice Boltzmann formalism; third, model 
potential or interaction, and free energies can be directly implemented into the lattice Boltzmann models 
[8]. It can be shown that the lattice Boltzmann equation recovers the near incompressible Navier-Stokes 
equations [9]. There is numerical evidence that the lattice Boltzmann method can indeed faithfully simulate 
the incompressible Navier-Stokes equations with high accuracy [10, 11, 12, 13]. 

Historically, the lattice Boltzmann equation evolves from its predecessor, the lattice gas automata [14, 
15]. Recently it has been shown that the lattice Boltzmann equation is a special finite difference form of 
the continuous Boltzmann equation [6, 7]. This result has set the mathematical foundation of the lattice 
Boltzmann equation in a rigorous footing [6, 7]. It also provides insights to relate the lattice Boltzmann 
equation to other existing gas kinetics based schemes [16, 17, 18]. In this paper, we discuss the mathematical 
connections between the lattice Boltzmann equation and the beam scheme [16]. 

This paper is organized as follows. In Sec. 2 we briefly discuss the BGK Boltzmann equation and its 
hydrodynamics. In Sec. 3, we describe the derivation of the lattice Boltzmann equation by discretizing 
the continuous BGK Boltzmann equation, and the derivation of hydrodynamic equations from the lattice 
Boltzmann equation via Chapman-Enskog procedure. We also analyze the validity and accuracy of lattice 
Boltzmann method. In Sec. 4 we describe the beam scheme and connect the beam scheme with the lattice 
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Boltzmann method. In Sec. 5 we compare the beam scheme and the lattice Boltzmann method. In Section 6 , 
we propose an improved beam scheme and conclude the paper. 

2. BGK Boltzmann Equation and its hydrodynamics. We begin with the Bhatnagar-Gross- 
Krook (BGK) [19] Boltzmann equation, which is a model kinetic equation widely studied [20, 21]: 

(2- 1 ) + = 

T 

where the single particle (mass) distribution function / = /(x, £, t) is a time-dependent function of particle 
coordinate x and velocity £, r is the relaxation time which characterizes typical collision processes, and g is 
the local Maxwellian equilibrium distribution function defined by 

(2-2) g(p, u, 6) = p (2tt 9)~ d ^ 2 exp [— (£ — u) 2 /20] , 

where D is the dimension of the velocity space £; p, it. and 9 = k B T/m are the mass density, macroscopic 
velocity, and normalized temperature per unit mass, respectively; k B , T, and m are the Boltzmann constant, 
temperature, and particle mass, respectively. The mass density p , velocity it, and the temperature 9 (or 
internal energy density) are the hydrodynamic moments of / or g: 

(2.3a) p = Jd£ f = (dig, 

(2.3b) pu = JdUf = JdUg, 

( 2 - 3c ) f P0 = = jdi\{i-u) 2 g. 

At equilibrium f ~ g, the Boltzmann equation ( 2 . 1 ) becomes 

(2-4) D t g = 0 , 

where 

( 2 - 5 ) D t =d t + i-V . 

The Euler equations can be easily derived from the above equation by evaluating its hydrodynamic moments: 


(2-6a) d t p + V- ( pu ) = 0 , 

( 2 - 6 b) d t (pu) + V(puu + p9) = 0 , 

pu 2 it) = 0 , 

where j = (D + 2)/2 is the ratio of specific heats of an ideal gas, ttit denotes a second rank tensor u t Uj, and 
u = it • it. It should be noted that 'y is related to the number of the degree of freedom of a particle, which 
is D in the case of mono-atomic ideal gases. The momentum equation ( 2 . 6 b) can also be rewritten as the 


following 


(2.7) 

pd t u + puVu = -VP, 

where 


( 2 . 8 ) 

P = p9 


( 2 . 6 c) 


D 


-p0 + 


ipit 2 )+ v -( 


7 pOu - 


2 



is the equation of state for ideal gas. 

The Chapraan-Enskog analysis [20, 21] gives the first order solution 

(2.9) f {x) = -rD t g. 

Therefore, 

(2.10) / = / (0) + / (1) = g - rD t g ■ 

With the above solution of /, the BGK equation, Eq. (2.1), becomes 

(2.11) D t g - TD\g = -1/ (1) • 

The moments of the above equation leads to the Navier-Stokes equations 

(2.12a) d t p + V-(pu) = 0, 

(2.12b) d t {pu) + V(/cmu + pO) = — ' , 

where l"^ 1 ) is the first order shear-stress tensor 

(2.13) n{J > = JdttitifM ^-rpeidiUj + djUi). 

It is obvious that the kinematic viscosity in Eq. (2.12b) is 


(2.14) 


v = rO . 


3. Lattice Boltzmann equation and its hydrodynamics. 

3.1. Lattice Boltzmann equation. The continuous BGK equation, Eq. (2.1), admits a formal integral 
solution given by [20]: 


f(x + £6t, £, t + 6 t ) = e" 4t/A /(*, t) + \e~ Sl/x J e l ' /x g{x + £f', £, t + t') dt ' . 


The lattice (BGK) Boltzmann equation can be derived by discretizing the above integral solution in both 
time and phase space [6, 7]. The obtained lattice Boltzmann equation is 


fa{x + e a $£, t -f ^t) — fa{ x i t) “ “~[/a(®* 0 fa 01 i 


where r = A /6 t , / a (x, t) = W a f(x , e a , i), /a eq) is the discretized equilibrium distribution function, and 
{e a } and {W a } are the discrete velocity set and associated weight coefficients, respectively. The discretized 
equilibrium distribution function fa^\ and both the discrete velocities {e Q } and their corresponding weight 
coefficients {W a } depend upon the particular lattice space chosen. For the sake of explicitness, we use the 
nine-bit lattice Boltzmann equation in two-dimensional space in the following discussion. In this case, we 
have W a = 2n0 exp(e^/20)w a , where 



( 4/9, 

a — 0 

(3.3) 

w a = 1/9, 

a = 1, 2, 3, 4 


[ 1/36, 

a = 5, 6, 7, 8 
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(3.4) 


e a = 


{ 


( 0 , 0 ), 

(cos <f> a , sin 4> a )c, 


a = 0 

a / 0, 


and <p a = (a - 1)tt/2 for a = 14, and (a - 5 )tt/2 + tt/ 4 for a = 5-8, and c = 6 x /6 t . (Here 0 = c 2 /3 for 
the nine-bit model has been substituted to obtain a uniform lattice structure [6, 7].) Then, the equilibrium 
distribution function of the 9-bit model is: 


(3.5) 




(eq) 


w a p 


3 (eg • u) 9(e a • u) 2 _ 3u 2 1 

c 2 2c 4 2 c 2 


The hydrodynamic moments of the lattice Boltzmann equation are given by 


( 3 ' 6a ) = 

a a 

(3.6b) = = 

a a 

< 3& > f - “) 2 /a - £ he. - u)V<«> ■ 

a a 

Note that the quadrature used in the above equations must be exact for these hydrodynamic moments in 
order to preserve the conservation laws [6, 7]. 

The algorithm for the lattice Boltzmann equation consists of two steps: collision and advection on a 
lattice space as prescribed by Eq. (3.2). The collision is accomplished as follows: first of all the hydrodynamic 
moments are computed at each lattice site {x} according to Eqs. (3.6) , the equilibrium /^ eq) can be calculated 
then according to Eq. (3.5). The distribution f a is updated on each site by using the relaxation scheme: 
/ a (x, t + St) — /a(®, t) — [/ a (x, £) — fa q ^(®> £)]/ T * After collision, f a advects to the next site (x + e a St) 
according to the velocity e Q , i.e., f a (x 4- e a 6 u t + 6 t ) = / a (x, t + 6 t ). It is obvious that the algorithm 
is simple, explicit, and intrinsically parallel. All the calculations are local and data communications are 
uniform to the nearest neighboring sites. 


3.2. Chapman- Enskog analysis. The hydrodynamics of the lattice Boltzmann equation can be de- 
rived via Chapman-Enskog analysis [20, 21] with the following expansion [9, 22]: 


(3.7a) 


(3.7b) 

where 


00 Xn 

/«(* + e ° s t, t+6 t ) = Yl o . 

n— 0 
oo 

/a = E e "^ ) - 

n=0 


( 3 - 8 ) D' t = (d t + e a -V) , 

and e = 6 t . The normal solution of the lattice Boltzmann equation, up to the first order in expansion 
parameter e (which is the Knudsen number), from Chapman-Enskog analysis is 

(3.9a) f(0) = /(eq) _ 

(3.9b) ,m = -rD'jm . 

With the above solution, we can, accordingly, derive the following governing equations for the lattice Boltz- 
mann equation through Chapman-Enskog procedure: 

(3.10a) D' t fW= 0, 

(3.10b) D’J^ -1(2 r - 1 )6 t D' t 2 fW = -- f £> . 

£ T 
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The hydrodynamic equations, the Euler and the Navier- Stokes equations, can be obtained by taking the 
moments of the above governing equations. 

3.3. Hydrodynamics of lattice Boltzmann equation. The moments of the zeroth order equation, 
Eq. (3.10a), in the discrete momentum space [defined by Eqs. (3.6)] lead to the Euler equations 


(3.11a) 

d t p + V-{pu) = 0, 

(3.11b) 

d t ( pu ) + V(/96> + puu) = 0 , 

(3.11c) 

d t (~p0 + \p u *) + V-(7p0u) = 0 , 


Note that Eq. (3.11c) differs from its counterpart, Eq. (2.6c), derived from the continuous BGK Boltzmann 
equation, Eq. (2.1). The energy flux \pu 2 u due to the advection of fluid is missing. The omission of the 
energy flux term | pu 2 u can only be justified if 

1 , 

-pn «7 pQu, 

where u = ||u||. Because 0 = <? s = c 2 /3 in the case of the nine-bit model, where c s is the sound speed, 
therefore 


y^= y/(D + 2j, 

where M = u/c s is the Mach number. From above analysis, it becomes clear that the lattice Boltzmann 
equation is only valid for low Mach number flow, or in the incompressible limit. This is consistent with the 
low velocity expansion made to obtain the equilibrium distribution function 

The moments of the first order equation, Eq. (3.10b), give the Navier-Stokes equations similar to 
Eqs. (2.12) with the viscosity given by 

(3.12) i/=l(2r-l)<5 t 0= l(2r-l)^, 

where 0 = c 2 / 3 and c = S x /S t . The Navier-Stokes equation from lattice BGK Boltzmann equation can 
be easily derived by noticing the similarity of Eq. (3.10b) and Eq. (2.11), provided that the quadrature to 
evaluate the moments 

(3.13) f dCC m e~ <2 = 5Z«?„C 

^ a 

is exact for m < 5, because rn — 1 from the first order moment, m = 2 from fa q \ and m — 2 from the 
term (e a ■ V) 2 /i eq) in Eq. (3.10b). To obtain the two-dimensional nine-bit lattice Boltzmann equation, the 
third order Gaussian quadrature is the optimal choice to the aforementioned goal of evaluating the necessary 
hydrodynamic moments exactly [6, 7]. 

In the above derivation of the Navier-Stokes equation, and the governing equation Eq. (3.10b) in par- 
ticular, one can immediately realized that the accuracy of the LBE method is of second order, as previously 
speculated [23], because all the second order terms in the Taylor expansion are included in Eq. (3.10b), and 
the truncation error is of third order. The term —S t 0 / 2 in the viscosity is the manifest of the inclusion of 
second order terms (in space). The simplicity of this proof for the second order accuracy of the lattice BGK 
Boltzmann method is due the simplicity of the collision operator in LBE method. 
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4. Beam Scheme. The beam scheme [16] is a finite volume, gas-kinetic based scheme to solve hy- 
drodynamic equations. In the beam scheme, hydrodynamic variables (mass density p. momentum pu, and 
temperature 0) are given at a particular time in each volume cell. The equilibrium distribution function 
constructed from the hydrodynamic variables can be approximated by a finite number of “beams,” or a 
distribution of finite number of discrete velocities. Consider in an one dimensional case in which we want to 
use three discrete velocities in the velocity space, then the equilibrium distribution g is approximated in 
coordinate with three Kronick delta functions: 


9x = P (2710) 1/2 exp{— (£ x - u x ) 2 / 20} 

(4- 1 ) ~ p[ao£(£x - tii) + ai<5(£x - u x 4- Au x ) + a 2 fi(£ x - u x - Att x )] . 

We can calculate the unknowns ( a o , oj, a 2 , and Au x ) from the following moment constraints: 


(4.2a) 

(4.2b) 

(4.2c) 

(4.2d) 

The equations for the unknowns are: 



(# = fcx (Cx - 

3 P0 2 = 


^x) 9x 

Ux) 4 9x ■ 


And the results are 
(4.3a) 

(4.3b) 

(4.3c) 


ao -b ai + a2 = 1 , 
(ai -a 2 )=0, 
(ai 4- < 22 ) Au^ = 6 , 
(a\ + 02 ) A u^. = 30 2 

2 

a ° ~ 3 ’ 

1 

ai = a 2 = - , 

A u x — \/3^ . 


Therefore, there are three “beams” or “particles” in the beam scheme with the velocity u x - y/3 0, u x , and 
u x + y/36. Note that the weight coefficients of these three particles, a 0 , and a 2 , are identical to those 
derived in the lattice Boltzmann equation by using Gaussian quadrature. Thus, in the situation of u x — 0, 
the beam scheme is similar to the lattice Boltzmann equation with the three discrete velocity of — c, 0, and 
c in one-dimensional case, where we have substitute 30 = c 2 for isothermal fluids. 

However, the difference between the lattice Boltzmann equation and the beam scheme outweighs the 
similarity between the two, for the reason that the lattice Boltzmann equation is a finite difference scheme, 
while the beam scheme is a finite volume one. In the beam scheme, the “particles” move in and out each 
cell according to the velocity of these “particles.” After this advection process, the hydrodynamic quantities 
(p, u, and 0) are obtained through an averaging process in each volume cell. The “particles” with different 
velocity are mixed first to compute the averaged hydrodynamic quantities in the cell, and redistributed 
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through the calculation illustrated previously. This mixing (or averaging) process inevitably introduces 
artificial dissipation. And this dissipation is implicit, just as for any other upwind finite volume scheme. 
Therefore, the transport coefficients, such as viscosity, cannot be explicitly derived in the beam scheme, and 
thus the beam scheme cannot solve the Navier-Stokes equations quantitatively. 

5. Lattice Boltzmann Method and Beam Scheme. We now compare the pros and cons of the 
lattice Boltzmann equation and the beam scheme. Theoretically, the lattice Boltzmann equation accurately 
approximates the incompressible Navier-Stokes equations [9]. The method is simple, explicit, and intrinsically 
parallel. The transport coefficients can be obtained explicitly, and therefore there is no numerical dissipation 
in the simulations by using the lattice Boltzmann method. In addition, the lattice Boltzmann method is a 
intrinsically multidimensional scheme. The disadvantages of the lattice Boltzmann method are the obvious 
consequences of the low Mach number expansion and the regular lattice structure. Because of the low 
Mach number expansion, the lattice Boltzmann method is limited to incompressible flows and therefore is 
not applicable to compressible flows and shocks. In addition, the regular lattice structure of the lattice 
Boltzmann method is a direct consequence of constant temperature, i.e 9 = c 2 / 3 =constant in the case of 
the nine-bit lattice Boltzmann model. This, we suspect, is partly the reason for the failures of the thermal 
lattice Boltzmann models. Furthermore, because the equilibrium distribution function is not a Maxwellian, 
the H ~ theorem for the continuous Boltzmann equation no longer holds for the lattice Boltzmann equation. 
That is to say, the is an attractor of the lattice Boltzmann equation, but it may not be the equilibrium 
in the sense of the //-theorem. 

In contrast to the low Mach number expansion used in the lattice Boltzmann equation, the Maxwellian 
equilibrium distribution is expanded around the averaged macroscopic velocity in each volume cell. The 
discrete velocity set depends on the averaged velocity and temperature within each cell, and therefore it 
varies from cell to cell. The beam scheme leads to correct hydrodynamic equations, including the energy 
equation. Consequently, the beam scheme can well capture thermal and compressible effects in flows. Thus 
it is suitable for simulations of high-speed (hypersonic) flows, and it is also more stable for high Reynolds 
number flows. The beam scheme is a finite volume, upwind shock capturing scheme. Its natural shortcomings, 
like any other such schemes, are that it has intrinsic and implicit numerical dissipations due to the mixing 
of particles in each volume cell, and the transport coefficients cannot be obtained explicitly. Therefore it 
cannot solve the Navier-Stokes equations quantitatively. 

6. Conclusion. As we have shown in this paper, while the lattice Boltzmann equation and the beam 
scheme shares the same philosophy in the discretization of velocity space (in one dimensional space) - all 
the conserved quantities are preserved exactly in the process of discretization, their distinctive difference lies 
in their equilibrium distribution function. The lattice Boltzmann equation expands the equilibrium at u = 0 
and uses a polynomial (of u) to approximate the Maxwellian, therefore the method is limited to apply only 
to the near incompressible Navier-Stokes equations. The beam scheme obtains particle beams around the 
average velocity of the Maxwellian distribution, thus avoiding the low Mach number expansion in the lattice 
Boltzmann method. Naturally, the beam scheme is suitable for shock capturing in the compressible flows. 
Moreover, the lattice Boltzmann equation evolves on a lattice structure, information advects exactly from 
one node to another thus there is no mixing process involved. In contrast, the particle beams in the beam 
scheme move from one volume cell to another, and the mixing among the beams occurs in the construction 
of local Maxwellian equilibrium. Because of the uncontrollable numerical dissipation in the beam scheme 
caused by the mixing, it is difficult to use the scheme to simulate hydrodynamics quantitatively. 

It is interesting to compare the lattice Boltzmann method and the beam scheme in multi- dimensional 
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space. In two-dimensional space, the beam scheme only uses five velocities [16], i.e., one central beam with 
the bulk velocity ( u x , u y ), and two side beams each in x- and y-directions, (u x ±V 36, u y ), and (u x . u y ±\/W). 
Based on the analysis of the lattice Boltzmann equation [9, 14, 15, 22], it is well understood that such a 
discrete velocity set inevitably introduces anisotropy into the hydrodynamic equations resulted from the 
scheme. To removed the anisotropy, one can use the nine velocity set derived in the lattice Boltzmann 
equation. That is, the diagonal velocities, (u x ± vW , u v ± \/30), must be included in the two-dimensional 
beam scheme. With this modification, the volume cell in two-dimensional space become an octagon, instead 
of a square. This is equivalent to use the product of two one- dimensional Maxwellian, each approximated 
by three “beams.” This is feasible because Maxwellian is factorizable in the Cartesian, or other, coordinate 
system. 
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